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Abstract. The 3-point vertices of QCD are examined at the symmetric subtraction point at 
one loop in the Landau gauge in the presence of the Gribov mass, 7. They are expanded in 
powers of 7^ up to dimension four in order to determine the order of the leading correction. 
As well as analysing the pure Gribov-Zwanziger Lagrangian, its extensions to include localizing 
ghost masses are also examined. For comparison a pure gluon mass term is also considered. 



1 



1 Introduction. 



In the late 1990's an interesting property of the running of an effective couphng constant derived 
from the triple gluon vertex in Quantum Chromodynamics (QCD) was revealed. In several 
articles, [H El EJ HI |5], using lattice gauge theory techniques and restricting to the Landau 
gauge the effective coupling constant appeared to deviate from the expected behaviour in an 
energy range intermediate between high and low. Moreover the deviation between expected and 
measured behaviour could be fitted by a power law correction. While such discrepancies are not 
ordinarily unexpected in essence because of our knowledge of the operator product expansion, 
the power correction was claimed to correspond to a dimension two operator rather than a 
dimension four one, [H EJ [3l HJ [5] . The latter is usually associated with the vacuum expectation 
value of the square of the gluon field strength, which is gauge invariant, and is termed the gluon 
condensate. However, it was proposed that the operator associated with the dimension two 
correction was ^^JJ^ where is the gauge potential. Such a gauge variant operator is not 
excluded as ultimately the running coupling constant is not a physical quantity. Subsequent to 
this there has been interest in trying to understand this dimension two operator both in the 
coupling constant situation and other quantities such as the gluon and Faddeev-Popov ghost 
propagators in the infrared. The condensation of a dimension two operator could be related to 
confinement in the sense that it generates an effective gluon mass. See, for example, [6] for a 
summary. Thus in the infrared the gluon propagator will freeze to a finite non-zero value at zero 
momentum. This is in accord with a long-standing result of Cornwall, [7J, who established that 
the frozen gluon propagator follows from the dynamical generation of an effective mass which is 
momentum dependent. In more recent years progress with Landau gauge lattice computations 
has produced data which actually appears to support a frozen gluon propagator. For instance, 
see the early activity in this respect PElIinilllliniliailllllSlIIH]. This should be qualified 
by the remark that ensuring that one is in a properly fixed gauge free from the complications 
of Gribov copies, [17j . is a non-trivial exercise. Indeed it is perhaps fair to comment that this 
has yet to be fully resolved. Also the zero momentum regime is numerically difficult to achieve. 
Aside from the gluon propagator the Faddeev-Popov ghost propagator has also been studied in 
the Landau gauge. It too can be used to study an effective running coupling constant at zero 
momentum via the ghost-gluon vertex which exploits properties derived from the Slavnov- Taylor 
identities, [18]. 

From the analytic point of view one can focus on the Landau gauge in the infrared. The 
seminal work in this area was by Gribov, [17], who highlighted the inability to fix the gauge 
globally due to the presence of Gribov copies. One property he established was that a less 
incorrectly fixed Landau gauge could be effected with a modification of the Yang-Mills action. 
This restricted the path integral to the first Gribov region and introduced a non-local operator 
into the Lagrangian. Although a semi-classical approach was used in [17] the non-locality was 
localized in a series of articles by Zwanziger, [ISlEOlEIlEaEallMlESlESlET], to produce 
a local renormalizable Lagrangian. This construction meant that one could compute in the 
Gribov context and study infrared behaviour. Already noted in [IT] the copies introduced a 
new mass called the Gribov mass, 7, which is not an independent parameter as it satisfies a 
gap equation. Its presence ensures that the Faddeev-Popov ghost propagator enhances and the 
gluon propagator is suppressed in that it vanishes at zero momentum. While this does not accord 
with recent lattice data, modifications of the original Gribov-Zwanziger Lagrangian, 128 1 129 ^ 130]. 
have been developed which do model the frozen gluon and non-enhanced Faddeev-Popov ghost 
propagator. Though such refinements are not unique. In |30j it was demonstrated that several 
different localizing ghost condensations can model a frozen gluon. However each has a different 
prediction for other quantities which are not yet or as widely measured on the lattice. While 
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there has been an intense amount of lattice activity on the propagators with as yet no resolution 
as to which refined solution is the leading candidate, the presence of an independent mass scale 
in the pure or refined Gribov-Zwanziger Lagrangians could be the source of the power corrections 
of the effective coupling constants of [U [2l O [H [5] . Therefore, it is the purpose of this article to 
explore this possibility and produce the one loop corrections to 3-point vertices of QCD using 
the Gribov-Zwanziger setup. We will do this for the specific momentum configuration of having 
the three external legs at non-zero momenta and with squared momenta all equal. This is known 
as the symmetric point and is a non-exceptional momentum configuration. Hence it will be free 
of infrared ambiguities which could plague the asymmetric subtraction point calculation which 
is an exceptional configuration. 

One advantage of this symmetric subtraction point is that by considering non-zero external 
momenta one will avoid having to go to the far infrared which is beset with potential gauge fixing 
issues. In essence the momentum range we have in mind for comparisons to numerical work is 
the intermediate one where the explorations of were centred. In other words one is in a next 
to high energy approximation where one can formally access power corrections. At much lower 
energies the explicit expressions for vertex amplitudes would be complicated functions of the 
masses and momenta. These would then have to be expanded in a power series to reveal power 
corrections. This intermediate energy range is important from the Feynman diagram point of 
view. The main reason for this is that we can expand the Feynman graphs in a power series in a 
mass scale. This avoids having to actually determine the explicit complicated functions and then 
expand them. Central to our calculations will be the use of the method of |31j which allows us to 
correctly power expand Feynman integrals. Given that the Gribov mass is the associated scale 
in the pure Gribov-Zwanziger case there is no a priori reason to exclude it as being a potential 
source for the power corrections advanced in [H El EJ IH EJ . For instance, one could in principle 
relate 7^ to the vacuum expectation value of precisely because 7^ appears in the gluon 

propagator, [T7|. Hence (^^^^) will be proportional to 7^^ on dimensional grounds but both will 
be as indistinguishable from the other as the dimension two correction measured on the lattice. 
Thus in our exploration of this problem the aim will not only be the determination of whether 
there are dimension two, four or higher corrections but also to deduce the magnitude and sign 
of the coefficient. In studying the triple gluon, ghost-gluon and quark-gluon vertices we will be 
able to devise a test which in principle could distinguish which of the refined solutions if any is 
favourable but crucially in a regime where deep infrared issues, such as Gribov copies, do not 
complicate a lattice study. However, as a control on the Gribov mass investigation we will also 
repeat the same analysis for the case where the gluon has an explicit mass, m. Ordinarily such 
a naive mass term breaks gauge invariance but one can have a non-local but gauge invariant 
gluon mass term. (See, for instance, [32l[33].) This operator reduces to ^m^^J^^ in the Landau 
gauge. Having such a control calculation to compare with the Gribov mass case is important as 
a gluon mass also mimicks a frozen gluon propagator. Thus that could actually be the source 
of the lattice computations. If that were the case then a symmetric vertex study could be used 
as a confirmatory test. 

The paper is organized as follows. In section 2 we discuss the general background to the 
Gribov-Zwanziger Lagrangian and its extensions which are needed for the computations we 
carry out. This includes the formalism we use to construct the three vertex functions at the 
symmetric point. Section 3 focuses on general aspects of the calculation including how the 
Feynman integrals are deduced in a power series expansion using the technique devised in [31j . 
The subsequent four sections are devoted to recording the results for the pure Gribov-Zwanziger 
case, what is termed the Q and Tt solutions and a pure gluon mass respectively. We discuss our 
test in the concluding section. Several appendices provide technical details which supplement 
the main discussion. 
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2 Background. 



We begin by recalling the main aspects of the Gribov-Zwanziger construction which are necessary 
for our computations. In [17j the Yang- Mills action was modified in order to take into account the 
ambiguities which arise in the gauge fixing procedure. With the restriction of the path integral 
to the first Gribov region the modification defines a boundary termed the Gribov horizon which 
manifests itself as a non-local dimension zero operator in the Lagrangian. It is related to the 
Faddeev-Popov operator, {d^^D^Y^ ^ in that the boundary corresponds to the surface defined by 
the first zeros of the operator. Thus the inverse operator will be infinite at the Gribov horizon 
and hence as long as there are no poles in the inverse one is within the Gribov region, [17]. 
Originally in [I7j, the action was treated in a semi-classical approximation in such a way that 
only the leading term of the Faddeev-Popov operator was used to define the boundary in the 
path integral. Subsequently, Zwanziger extended the analysis to all orders to produce the inverse 
Faddeev-Popov operator, [21]. The resulting Lagrangian is, [26j . 

^Gribov ^ ^QCD ^ ^;eac^e6d^a(^)|^_^y'^fcM(^) _ (2.1) 

where g is the coupling constant, d is the spacetime dimension, A'^ is the dimension of the 
adjoint representation of the colour group whose structure constants are f"-^'^. Here 

^QCD ^ _ 1(7^ (^a^i^ _ _}_^Q^lJ^a.2 _ ^qi^j:, + i^^^^^^ (2.2) 
4 ^ 2a 

is the usual QCD Lagrangian valid at high energy where c" is the Faddeev-Popov ghost and 
ip'^^ are massless quarks. Although we have included the gauge parameter a associated with the 
linear covariant gauge fixing we will perform all our calculations in the Landau gauge which is 
a = 0. In (j2.ip the mass parameter 7 is known as the Gribov mass and is not an independent 
parameter as it satifies a gap equation derived from the horizon, ()2.3p . Only when 7 satisfies 
the gap equation is one actually in the gauge theory jTTl [211 ES] • The field independent term of 
([2.ip ensures the non-triviality of the condition since, [26], 



feacjeMl^a,^l^\ \ (2.3) 




9^ 



The presence of the non-local operator implies that the gluon propagator is modified from that 
which is used at high energy. Specifically, \n\ . 



{Al{p)K{-p)) = - J^^Iyfc^f^^^iP) (2-4) 



which behaves as — ^^^^^5°'^ at large momenta where 



PAp) = - (2.5) 
At low momenta the propagator tends to zero which is the gluon suppression, |17] . 



While (|2.ip extends the Yang-Mills action to incorporate the copy issue, from a practical 
point of view the non-local term means that one cannot use it for explicit calculations. To 
circumvent this Zwanziger managed to localize the non-locality in several articles, [20 [ I21 [ [26]. 
to produce a local Lagrangian. Its renormalizability was established in [26 [ I34 [ [55]. In order to 
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achieve the locahzation locahzing ghost fields were introduced, [26], so that (12. ID is replaced by 



^GZ ^ ^QCD ^ ^pab^^Qu^jj^p^^ab ^ ^_pab^.QU^JJ^^^ 



1 



\ab 



Z7 



2 jabc j^a ^i^bc 



dNAl^ 



2g^ 



(2.6) 



Here /9^^ and are real fields and we have chosen this version over the complex fields of the 
earlier localization, [2H I26j . Accompanying these bosonic fields are the other localizing 
ghosts, ujI^ and (D^*, which are Grassmann. They are required to ensure that the ultraviolet 
structure of the theory such as asymptotic freedom is not upset nor renormalizability lost by the 
sole presence of the bosonic ghosts. The other main feature is that the dimension zero operator 
of ()2.ip is translated into a dimension two mass-like term which mixes A'^ and In effect this 
term corresponds to the non-local term which can be clearly seen by recalling the relation 



4" 



be 



(2.7) 



which is deduced from the ^^'^ equation of motion. The original horizon condition (|2.3p becomes 



bcr 



21 



9' 



in (j2.6p . With the extra fields there is an extended set of propagators which are, [26], 



{C\p)C'i-p)) 



[(p2)2 + CA7l 
^ac^bd 



jabe jcde^i 



PAp) > {A't{p)pl%-p)) 



{pf{p)Pu{-p)) = {ooI\pW{-p)) 



p2[{p2)2 + CAl^] 
^ac^bd 



pz 







(2.9) 



The suppressed gluon propagator emerges again. Despite the presence of the mixed propagator 
it is possible to compute one and two loop corrections to various quantities such as the gap 
equation for 7, [iTl I2HI [36] . 

While this was the standard Lagrangian used to incorporate Gribov copies it did not cover 
the observed behaviour on the lattice since the gluon propagator does not freeze to zero but to 
a finite non-zero value. To account for this in the Gribov-Zwanziger context various extensions 
to (|2.6p have been considered in |28[ l3Uj . Each generalization can be summarized in the 
addition of a mass operator for the localizing ghost sector. The most general such operator is, 
[30], 



O 



where 



^2 .ac^bd ^ ^2^pcejbde ^ jabe jcde ^ ^2^^abcd ^ [4^ab^cd ^ ^^ad^bc 

Ca Na 



O 



abed 



pOhpCd ^ -^abpCd _ -pab^cd ^ ^ab^cd 



^abcd 

(2.10) 
(2.11) 



Each colour structure is tagged with a mass, px, where we use the same labels and conventions 
as were used in |30] . The tensor d'^'^'^ is totally symmetric and is defined by, |37] . 



j^abcd 



A-^A ^A-^A 



(2.12) 
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Overall the additional operator is BRST invariant and satisfies a Slavnov- Taylor identity so 
that the operator renormalization constant is related to the wave function renormalization 
constant, [281 [291 [30] . With this additional term the propagators (12. 9p become significantly more 
complicated. These were analysed at length in [30] where the full set for SU{3) were recorded 
explicitly. In general and for certain specific cases the extra mass parameters nx lead to gluon 
propagators which freeze to a non-zero finite value. In [29] one of these cases was examined 
which in the notation of (I2.10p was termed the Q solution. However, in pOl it was noted that 
this solution was not unique and that another specific single mass parameter solution could be 
equally viable. This was the TZ case. It was argued in |30j that this was a more natural solution 



than the Q case since a non-zero fi-ji would correspond to a condensation of the operator ^'^^^^-^t^ 
consistent with the structure of the pure Gribov-Zwanziger propagators, |3(3]. To see the frozen 
gluon propagator at the outset the propagators for these two specific solutions are 



and 



{A;{p)Ai{-p))Q 



{Pfip)pt'i-P))Q 

{A'^,{p)At{-p))n 
{Al{p)e%-P))n 

{e^mi\-p))n 







pAp) ^ {C{p)pf{-P)) 



[{p^Y + pIp' + Cai^]' 



Q 



+ pq] 



jabe jcde^A 



[p^ + plEp^f + pIp^ + Cai^ 



■pAp) 



(2.13) 



[{p2)2+^2^p2 + CAl^, 



■p. 



fJ,U 



{p) , {Al{p)pl\-p)) 



ifabc^2 



[{p'y + p'nP' + CAi'] 



pAp) 



{e^{p)pl\-p))n 



ac xbd 



jabef 



+ Cai' 



P 

jabe jcde 



CAP^ip^r+plp^ + CAY 



-.pAp) 



Pn 



{pt{p)p'if{-p)) 



n 



Cap^V + p\] 

{co^\p)Cot''{-p))n 



2^Lf,^{p) 



aczbd 



-r],iu + 



fabefcde^2^ 

Cap'^Ip'^ + 



V^u ■ (2.14) 



It is these propagators which we will use as part of our study of the power corrections to the 
3-point vertices. While the gluon sectors of each are formally equivalent the key differences are 
in the localizing ghost propagators. In the Q case there are no massless factors in any of the 
propagators whereas there are massless modes in the TZ solution. In |30j it was noted that this 
led to different infrared properties of the localizing ghost propagators after the gap equation for 
7 is satisfied. Such differing behaviour can be used to distinguish from these solutions if lattice 
data was available for the localizing ghost propagators. 

We now turn to the formalism relating to specific Green's functions we will compute in a 
power series expansion. These are 



Al{p)A',{q)Al{r: 



c'^{p)c\q)Al{r] 
i^\pW{q)Al{r] 



^J:^'^Hp,q,^^^^^) 



(2.15) 
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where 

r = — p — q (2-16) 

and each Green's function corresponds respectively to the triple gluon, ghost-gluon and quark- 
gluon vertices. The independent external momenta flowing in through two of the external legs 
are p and q and neither are nullified so that we are in a non-exceptional momentum configuration 
which does not suffer from infrared issues. All squared momenta are held at the same value 



p 



1^ 



which implies 



1 



pq 



(2.17) 



(2.18) 



where fi is the mass scale introduced to ensure that the coupling constant remains dimensionless 
in d-dimensional spacetime as we will be using dimensional regularization throughout. Our 
regularizing parameter will be e where d = 4 — 2e. As each Green's function carries colour and 
Lorentz indices we have to decompose them into scalar amplitudes. For the former we have done 
this in (|2.15p . This is relatively straightforward since to the loop order we are working there is 
only one colour tensor for each vertex function which is evident from the explicit computations. 
For the Lorentz sector we have to introduce a set of basis tensors which are built from r/^j^ and 
p^ and q^. Thus 



(p,9,7^/^x) 



14 



fc=l 

2 

k=l 



[k) 



(k) 



(p,g,7^/xl) 



fc=i 



(2.19) 



which defines the scalar amplitudes for each vertex function. The explicit tensors for each case 
are given in Appendix A and we use the same set as was used in |38j. We note that away 
from the symmetric point restriction the basis will involve more tensors. Also the basis is not 
unique and there are other choices. It turns out that for the triple gluon vertex from the explicit 
computations of [38j we have checked that to two loops one can write the vertex function more 
compactly in terms of three tensors. One of these, for instance, corresponds to the Feynman 
rule of the vertex itself when r = — p — q. The other two do not involve rj^i, and their explicit 
forms are given in Appendix A. Therefore, here we choose to work in this more compact basis 
and replace the first equation of (j2.15p by 



3 



/ ^ (fc) fj,ua 



{p, q) s 



{k) 



(2.20) 



To distinguish this basis from the previous set the amplitudes will have a tilde. In order to 
calculate each scalar amplitude we use the same projection approach as [38]. Briefly each 
Green's function is multiplied by a linear combination of the basis tensors in d-dimensions. The 
coefficients are found by first constructing the matrix of products of all the basis tensors. Then 
the inverse of this matrix gives the linear combination. Denoting this inverse matrix by M we 
have 



eabc 



f''-1^1^f[p,q,r 

r'^^TkUp^q,^', 



\k) 



A) 



A) 



M 



kl 



^^^^■"'{p,q) {A';^{p)Ai{q)A^^{r] 



Ml? 



{rl^^^iP^q) {c''{p)c\q)Al{r: 
(^3f''(p,g) (i^\pW{q)AUr: 



(2.21) 
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where k and / are the matrix labels. The explicit matrices are given in Appendix A. For the 
quark sector we have the additional issue of spinor indices to account for in the amplitude 
decomposition. So as well as building the tensor basis from r]^^, and (7^ one has also to 
include 7^. This implies that products of 7-matrices can arise in the tensor basis and since we 
will be working in d-dimensions it is natural to use the generalized basis of 7-matrices which 
spans spinor space in d-dimensions, |39l SOI HI]. These are denoted by r|'^j"^" where n is a 
positive integer and defined by 

pjn...Mn ^ ...7/^"] (2.22) 

where 1/n! is included in the antisymmetrization. The algebra of these has been studied at 
length in [HI US]. Though for constructing the matrix A^^J^^ a useful property is, [32| 

tr(rf^y'^'"r^^y'^") oc (2.23) 

which partitions the matrix where /Mi---Mmi'i...i^n jg i^j^g y^nit matrix on the generalized spinor 
space. 



3 Calculational method. 

We now detail the overall method we have used to determine the power corrections to the one 
loop vertices at the symmetric point in each of the various Lagrangians we are interested in. 
We have followed a general approach which allows us to construct routines for all possible mass 
configurations. For instance, examining the propagators all possible one loop 3-point functions 
will have at most three different non-zero mass scales. However, since we are at a symmetric 
point and there are massless poles in the propagator sets, there are actually only seven distinct 
basic mass distribution configuration we have to consider. If we define 



2 2 2 



{2ttY [A;2 + mfl [{k - + m|] [{k + qf + ml 



(3.1) 



then the basic master integral structures from the point of view of mass distribution across the 
propagators are 1(0,0,0), /(mf,0, 0), /(mf , , 0), /(mf,m2,0), /(mf , , ), I(mf , , 7712) 
and /(mf, ml, 771.3). Here none of the rrii are equal. We note that (j3.ip is the basic structure 
since any of the Gribov or Stingl type propagators can always be written as the product of 
two canonical propagators and then partial fractioned. The fact that we are at a symmetric 
subtraction point means that I{mi,m2,Tn1) is totally symmetric in its arguments. We have not 
included powers of the propagators in the definition of I{m1,m2,ml) as we use the standard 
approach of breaking up all the contributing Feynman graphs into scalar integrals. These are 
determined after the general projection method. Then numerator scalar products are rewritten 
in terms of the propagator factors. For one loop 3-point functions at the symmetric point there 
are no irreducible numerators. To proceed further we have to write these integrals, which may 
have propagators to a negative power or positive power greater than unity, in terms of the basic 
master integrals for each of the seven mass distributions. For the massless case we did this in [38j 
at one and two loops. For the massive cases we extend that approach which used the Laporta 
algorithm [U]. That method, [H], uses integration by parts to establish algebraic relations 
between all the integrals which arise. As there is an overredundancy in the relations, one can 
relate all integrals down to a basic set known as masters. The values for these are determined 
by other methods and hence the whole calculation is complete. In terms of tools we have used 
the Reduze implementation, [35], which uses the GiNaC symbolic manipulation system, |46j . 
and is written in G-|— |-. The Reduze package creates a database for each topology and mass 
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distribution and one extracts the relations which are needed for the specific computation. Aside 
from the massless one used in [38] we have constructed five other one loop databases. That for 
/(mf, 771.2, ^-i) is needed as no integral of this mass distribution arises for propagator powers 
other than unity. Once the databases are determined the required integral relations are written 
in Form, |i47j . and included in the automatic computation routine. To generate the Feynman 
graphs for each of the three vertices we use Qgraf, [l8]. For the Gribov-Zwanziger case there 
are 30 one loop graphs and for the other two vertices there are 3 one loop graphs in both cases. 

Having broken down the computation of the amplitudes of each vertex into the basic master 
integrals we have to substitute the explicit values for each. In general terms at this stage there 
are three basic configurations which correspond to having one, two or three propagators. By 
this we mean that in the latter two cases any two or three of the propagators present in (j3.ip . 
For the one and two propagator cases for different masses these integrals are known exactly 
and their power series expansion can then be substituted after carrying out a Taylor expansion. 
For I {nil, 1712, m'^) itself the explicit result is not known for any combination of non-zero 
except when all three masses are zero. Therefore, since we are only interested in the power 
series expansion we follow the method of [31]. This method allows one to expand Feynman 
integrals in powers of rr? j [i^ where m is a generic mass scale deriving from a propagator in the 
original integral and [i is our common scale here for the squared momentum of the external legs. 
While [31j detailed at length the expansion of the two loop self-energy topology as an example 
the method is general and we use the general formalism that was provided there. If in general 
we denote by J one of our master integrals /(mf , m^, m-3) with at least one ^ then the 
asymptotic expansion is, [3T] . 

Jr ~ ^ Jvi\°T\m,Y\q,}'K ■ (3.2) 

A 

We use similar notation as [3_lJ and note that F is the original Feynman diagram, A are subgraphs 
which arise in the asymptotic expansion and rrii formally represent the masses on each propagator 
of F. (In [31j 7 was used for the subgraphs but we use A here to avoid confusion with the Gribov 
mass which is a parameter in the actual expansion.) These are necessary to counteract the 
infrared infinities which arise if one naively expands the original integral in powers of 1//U^. 
That is always the first term in the expansion and in that case the subgraph A is the unit graph. 
The other non-unit graphs in the sum are constructed from all possible routings of the external 
momenta around the graph. In the two loop example detailed in ^31] there was only one such 
momenta unlike the two here. However, it is straightforward to see that there are three such 
graphs since with two momenta there are three ways to route the external momenta. In other 
words in each of the three graphs one of the three propagators of (13. ip will have no external 
momenta. For each of these three cases the subgraph A is expanded in powers of the masses 
rrii and the momenta qi external to that subgraph itself. The terms of this Taylor expansion, 
denoted by 7{mi};{gi}'^Ai are then substituted in the reduced diagram J^jx before performing 
the integration over the loop momentum, [31j. For the expansion the leading term will involve 
massless 3-point integrals which can be reduced to the one loop massless master of jl9j using the 
Reduze database we have already constructed. The remaining three terms reduce to one loop 
massive tadpoles which are readily evaluated. In Appendix C we have given various examples 
of the expansion of master integrals which were used within our calculations. 

4 Power corrections in the Gribov-Zwanziger Lagrangian. 

We now turn to the mundane task of recording the explicit results for each of the three vertices 
for various versions of the Gribov-Zwanziger Lagrangian, its extensions and the gluon mass 
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cas^. In this section we focus on the pure Gribov-Zwanziger Lagrangian, (I2.6p . where there 
are no mass terms deriving from locahzing ghost field masses. Throughout this and subsequent 
sections we give the symmetric vertices to and including the term which corresponds to the 
dimension four correction. While we could in principle include higher order corrections, there 
does not appear to be any practical reason to do this at present since such terms would be 
difficult to extract numerically from the lattice. First, we record that the triple gluon vertex 
structure for (12.611 is 



71^^,9, 7^0) 



(1) 



+ 



13 77r 



7 



6 36 24 

m^{p,q, 0,0) 



2 



Ca^ 



+ 



37r Cfj^ 



32 //2 

m^'{p,q, 0,0) 



+ 



289 787 , 

\ In 

144 192 



+ o 



+ 



32 /i2 



+ 



25 
18 



25 2 2033 799, 

— TT^ \ In 

27 576 192 



a + 0(c? 



a + 0{a^ 



(4.1) 



Interestingly there is no dimension two correction for channel 1 which is the channel correspond- 
ing to the original triple gluon vertex Feynman rule. The other amplitudes have a dimension 
two correction and moreover, that term is the same in both cases. The explicit values of the 
7^ = amplitudes for this and the other cases are given in Appendix B. For the ghost-gluon 
and quark-gluon vertices we have 



m^{p,q,7^0) = J:ll^{p,q,0,0) 

32 

^^^{p,q,l',0) = J:T2^{p,q,0,0) 
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For the latter the equivalences between various amplitudes merely reflects the underlying left- 
right symmetry of the vertex in the choice of basis tensors we have used. We have not imposed 

this symmetry within the computation but instead it has emerged naturally and is regarded as 
a minor internal check on the programming. Unlike the triple gluon case there is a dimension 
two correction for all channels of both these vertices. 
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5 Q solution. 



Next we turn to what is termed the Q solution in the notation used in [29J. In this and 
subsequent sections the results are more involved due to the presence of an additional mass 
scale corresponding to the localizing ghost field mass of the particular solution. Therefore, we 
will have various combinations of 7 and fix- To compactify notation we will define 
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where I corresponds to the particular solution of interest. For the triple gluon vertex we have 
the expressions 
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where like the pure Gribov-Zwanziger case the structure of dimension two corrections are the 
same. This is somewhat unexpected as the presence of the extra mass does not induce a lower 
order correction. Therefore, if one were using the absence of a dimension two correction in 
lattice data in channel 1, if for example that were the case, then one could not conclude that 
the pure Gribov-Zwanziger case is an explanation since Q has the same qualitative feature. We 
note that here and subsequently we have used 0{iJ,j/^^) to indicate we are dropping dimension 
six terms and have used //j within the order symbol to indicate a generic dimension six mass 
scale. In practice this could be a dependent on 7 too but it would complicate the notation. For 
the other two vertices we have 
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Again these expressions all have dimension two corrections albeit complicated. However, they 
have the same qualitative structure as the pure Gribov-Zwanziger case. In the above expressions 
for all the amplitudes we have checked that the results of the previous section are reproduced 
as /Xg — > 0. This also applies to results in subsequent sections in the appropriate /Ltj — )• limit. 
We have also checked that in the 7^ — limit there are no corrections at all as there should be 
since then there is no horizon condition. 



6 71 solution. 

Structurally, the expressions for the TZ are very similar to those for Q. For instance, in many 
cases the differences are only in the numerical coefHcients. Though the mass is fi-ji rather than 
fiQ of course. Therefore, we will make minimal comment on these parallel results. For the triple 
gluon vertex we have. 
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Again there is no dimension two correction for channel 1. For the other two vertices it transpires 
that aside from mapping the masses fiQ f-)- uti the expressions for the ghost-gluon and quark- 
gluon vertices are formahy the same at one loop. This is not unexpected if one considers the 
contributing Feynman diagrams. The propagators for Q and TZ differ only in the localizing 
ghost sector. Since the gluon propagator is the only propagator from the gauge sector which 
contributes to the ghost-gluon and quark-gluon vertices at one loop and is formally the same 
for both solutions then the respective vertex functions have to be the same. 



7 Explicit gluon mass. 

So far we have concentrated on Gribov and Stingl type propagators for the gluon sector of Yang- 
Mills based on the pure and refined Gribov-Zwanziger Lagrangians. However, an alternative 
scenario is that a frozen gluon propagator in the infrared could be as a result of a pure gluon 
mass. While such a term ordinarily breaks gauge invariance, it is possible to construct a gauge 
invariant mass operator for the gluon. However, such an operator has to be non-local, [32> ,33] . 
and could be associated with an analogous ghost operator which is non-local. The origin for the 
latter observation is that it is possible to have a BRST invariant gluon mass term which was 
considered in [50J. There in order to make that local gluon mass operator BRST invariant one 
had to include a ghost mass term which was dependent on the gauge parameter. In either case 
when one restricts both to the Landau gauge the non-locality disappears in the first case and 
in the second only one of the two terms contributing to the BRST operator survives. In both 
cases in the Landau gauge 

Oai = ^A^A-^ (7.1) 

emerges as the mass term for the gluon. We use the same notation as |30j for consistency and 
will regard the presence of (|7.1|) in a Lagrangian as our control calculation. With (|7.ip we can 
derive the propagator for j4^. However, in order to compare with the previous two cases we will 
also add (|7.ip into the Gribov-Zwanziger Lagrangian. Thus our propagators are 
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Clearly in the limit where 7 — )• the first term reduces to what one expects for a gluon propagator 



with a naive mass term. Although the , ^'^ and u^f propagators remain their contribution 
cancels within any calculation as if they were not present in the first case in this limit. Though 
for any results we present we have checked this directly by using the usual Yang- Mills Lagrangian 
without including a Gribov mass term. 

The procedure to compute the vertex functions for this case is precisely the same as the 
previous two sections. However, we first give the results for a pure mass term, fj,x, in the 
absence of 7. For the triple gluon vertex we have 
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In contrast to the pure Gribov-Zwanziger case and the Q and TZ extensions there is a dimension 
two correction in channel 1 as well as the other two channels. This is a significant departure 
from the three Gribov scenarios and will form part of our test. For the other two vertices we 
have 
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Both these cases are completely parallel to earlier sections. 

Considering the situation where there is a Gribov mass as well as a gluon mass term we have 
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for the triple gluon vertex. The dimension two contribution in channel 1 remains as the leading 
correction but interestingly it is independent of 7. This is not the case for the other two channels 
in that the dimension two term derives from a combination of the two basic scales jjix and 7. 
So it would appear that at one loop there is a clear way of distinguishing between the masses 
in the problem. If a lattice computation found a dimension two contribution in channel 1 then 
that would be a clear indication of an explicit gluon mass term. Though it would not determine 
whether or not there was also a Gribov mass present too. That could be deduced from the 
details of other channels and vertices. For the ghost-gluon and quark-gluon vertices we have 
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These are considerably more complicated than their Q and TZ counterparts. However, the key 
point in this respect both here and in previous results sections is that it is the actual mass 
scale, no matter how complicated it is in each underlying theory, which is ultimately what the 
lattice will observe. We note that as a check for each of the three vertices the pure gluon mass 
expressions are obtained in the 7^ — >• limit in each case. 



8 Discussion. 



We conclude by recalling that we have analysed the one loop corrections to the 3-point vertices 
of QCD at the symmetric subtraction point in the Landau gauge using the Gribov-Zwanziger 
Lagrangian and its extensions which incorporate the Gribov problem. The main motivation was 
to examine the power corrections to the amplitudes in order to see whether the leading correction 
was dimension two or four. While [1] focused on an effective coupling constant derived from 2 and 
3-point functions and examined the deviations from expected behaviour, we have concentrated 
on just the vertices themselves. This is because the effective coupling constant definition involved 
the behaviour of the propagator form factors which have been shown in the Gribov-Zwanziger 
case to have dimension two corrections, [51] • These would therefore dominate in an effective 
coupling and the leading vertex correction behaviour would not be distinguishable. This is 
important since we have shown that in certain amplitudes the leading correction is dimension 
four and not dimension two. Indeed for any asymmetric momentum configuration for the three 
vertex functions the leading correction is always dimension two. The particular cases where 
dimension four is leading is in the channel of the triple gluon vertex which corresponds to 
the Feynman rule of the vertex itself. This is the case not only for the pure Gribov-Zwanziger 
Lagrangian but also for its extension to either what has been termed the Q and IZ solutions. For 
a simple gluon mass in the absence of the Gribov mass the triple gluon vertex has a dimension 
two correction. Since we are at a symmetric subtraction point where the common scale of the 



24 



external legs is not small then examining the deviation from expected perturbative behaviour 
could provide an important test of the underlying mechanism. For instance, if in measuring the 
triple gluon vertex channel 1 a power behaviour deviation of dimension two is found then that 
would rule out a pure Gribov-Zwanziger or Q and Tt explanation. It would not necessarily imply 
that a pure gluon mass is the underlying reason. This is because there are more complicated 
extensions of (|2.6p , [30j , not considered here where those propagators could mimic the dimension 
two correction. We have not introduced these here as we believe of the full set one of Q or 7^ 
is naturally favourable as discussed in [30]. By contrast if a triple gluon vertex measurement 
indicated that the leading correction was clearly dimension four then that would suggest that 
the Gribov mass is playing a role. To decipher whether it is the pure Gribov-Zwanziger case 
would require examination of the strength of the relative corrections and also data on the other 
vertices. Though this clearly is at a level of fine detail. 

It is worth commenting on the current status of lattice measurements of the three vertices 
at the symmetric point. First, we note that there is only a small amount of data for this point 
compared with the asymmetric point. For the triple gluon vertex the original indication of a 
dimension two correction in the effective coupling constant was carried out in [Ij. However, in 
keeping with other analyses of this vertex and the other two it transpires that signals suffer from 
more noise than the corresponding asymmetric vertex. This is despite the fact that the latter 
configuration requires a zero momentum limit. Therefore, at this stage one can not yet make any 
meaningful contact with data on this vertex function in order to see deviations from expected 
perturbative behaviour aside from the original observations of [l]. Moreover, we note that in 
[1] the vertex function for the triple gluon vertex was decomposed into only two independent 
tensors and not three as we have done here. Appendix A provides more details on this point. 
So even if a comparison with data could be made with our power corrections it is not clear 
whether this would be meaningful given that the bases are different. Other studies of this vertex 
include [521 153j . In the former the four dimensional data, while noisy, show a general decrease 
towards zero momenta which is reinforced in the latter article. Thus again in these cases a 
direct comparison with power corrections is not currently viable. For the ghost-gluon vertex 
the data of [53] does not suffer from as much noise as the triple gluon vertex case. There the 
main observation is that the ghost-gluon vertex is effectively equivalent to the tree value with 
a small maximum about IGeV, [53j. For the quark-gluon vertex there is the additional issue 
of quenched versus dynamical data. The analysis of [54] was in the quenched approximation 
and at the symmetric point the data indicate a smoothly rising vertex function. Though the 
data is noisier than the asymmetric results presented there too and not sufficient in order to 
perform a comparison with power corrections. Despite this one hope is that with the advances 
in lattice technology in recent years the focus could return to all of the 3-point functions now 
that there seems to be a consensus on the zero momentum behaviour of the gluon and ghost 
2-point functions. 

Finally, we should qualify our overall remarks by noting that we have performed the analysis 
at one loop. One could regard this as a next to high energy expansion. However, there is 
no a priori reason why the leading dimension four correction of channel 1 of the triple gluon 
vertex should persist beyond one loop. The leading two loop correction to this channel could 
be dimension two. So a more careful test could be that if the leading correction is dimension 
two but relatively weak compared to the other leading one loop dimension two corrections in 
other channels then that could be evidence for a non-pure gluon mass explanation. To carry 
out a two loop extension to the expansion is in principle possible but is beyond the scope of 
the present paper. The technical calculational tools are clearly available. Though one would 
have a substantial number of Reduze databases to build for the two main topologies that occur 
at two loops, [38], to cover all the non-zero mass distribution possibilities. While a one loop 
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observation is by no means a proof it is intriguing that of the three 3-point vertices of QCD it 
is actually the fully symmetric one in terms of fields when examined in this power expansion 
specifically at the symmetric subtraction point which should have dimension four as the leading 
correction. 
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A Tensor basis. 

In order to assist with the interpretation of the results in this appendix we record the explicit 
forms of the tensors in the bases for each of the vertices. For the triple gluon vertex the tensors 
of the original of [38] are 

^fioW^^'^) = ^'i'^P'^P- ' ^fnw^^''^) = ^Pi^luQo , Tfu)^^AP^(l) = ^Qf^P-^Q^ 

'^tli)^.v^P^1) = JoQt^Q'^P^ , T^fu)^,uaiP^Q) = \my(l<r (A.l) 

where the first six correspond to the terms of the original vertex in the Lagrangian. How- 
ever, from explicit calculations, [38], it transpires that the basis can be compactified into three 
independent combinations which we define to be 

'^ll)^lua (P^ l) = ^l^'^P'^ ~ '^l^'^l'^ ~ "^"ni^aPu - Vat^q,y + VuaPt^ + '^Viy^qfi 

'^{iftj.ua (P^ l) = i'^Pf^P'^P'^ + PfiQiyPrr - PfiQuqa + 2qf,PuPa - 2qf,puqa - 2g^g^g^] ^ 
'^{sffiuAP^ 9) = [PiJ.Puqa - qtiPuPa + qf^Pi^qa - qtiquPa] ^ • (A.2) 

The first again is the triple gluon vertex when one sets r = — p — q. In the analysis of [1] only 
two tensors are defined. One corresponds to the first and the second is a linear combination of 
the other two. More specifically the second tensor of yy is proportional to 

-PfZ'^AP^^) - = - l^.il-rUr-pUp-q).. (A.3) 



Though the explicit two loop calculations of [38] indicate that this choice of basis is too limited. 
In order to project out the amplitudes the projection matrix for ()A.2p is 



P'i 27{d - 2) 



/ 3 -6 

16(d-2) 8(d-2) I (A.4) 
^ -6 8{d-2) 4(4(i-5) 



where the subscript denotes the symmetric point. The corresponding 14 x 14 matrix for 
'^f^^i/a ^P^ ^) ^^^^ ^® those for the quark and ghost vertices used here were given in [38] 
but, for completeness here, the tensors of the latter two respective bases are 

rl^f^iP^Q) = P. , ^(2vb,'?) = 'Z. (A.5) 
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and 



t'^^Ap^^) = ^ > t'^^Ap^^) = ^ > T'^liP^^) = ;^r(3).,,. (a.6) 

We use the convention that when an external momentum is contracted with a Lorentz index 
then the dummy index is replaced by that momentum. 



B Leading order amplitudes. 



In this appendix for completeness we record the explicit values for the various one loop am- 
plitudes of each vertex at the symmetric point in the absence of the Gribov mass. They were 
computed originally in [55] . We have 



E^ff(p,(7,0,0) 



S^|f(p,(7,0,0) 



16 /I 



—Tp' - - 2 vr 



^+ .27 
+ Oia') 

16 128 2 64 ,, /I 

\ vr^ ih' - 

27 243 81 \3 

+ 0(a2) 

64 2 28 32 ,, /I 

TT^ th' - 

243 27 81^ V3 
+ 0(a2) 



32 



81 



TpNf + 



3 23 2 23 ,, /I 

- H vr^ ^' - 

8 162 108^ V3 



97 67 2 67 ,, / I 

vr^ H - 

108 243 162^ V3 



67 56 2 28 ,, /I 

vr-^ H tb' - 

54 243 81^ \3 



Ca 



(B.l) 



for the triple gluon vertex. We note that we have checked that these expressions agree with those 
derived in [38] after converting the amplitudes of [38j to the basis used here. The amplitudes 
for the ghost-gluon vertex are 



S(Jf(p, 9,0,0) = -1 + 
S(2J(P,'?,0,0) = 
and those for the quark-gluon vertex are 



5 2 5 fl 
:7r — —yj 



72^ V3 



108 
1 



72 



5 2 1 



CAa + 0{a^ 



108 



Cau + Oia^ 



(B.2) 



m^{p,qAO) 



1 + 



13 13 2 13 fl 

\ TT^ ■0 - 

4 54 36 \3 



Ca + 



9 V3/ 27 



Cf 



+ 0{a' 



7 5^5 ,,fl 

- H vr^ V - 

3 27 16 V3 



Ca + 



9 \3j 27 3 



Cf 



a + 0{a^ 



^1l)iP^QM 

2 K..>f^ 



(4) 

- H vr^ V , 

3 27 9^ V3 

27 18^ \3 



Ca - -Cf 



a + 0(a2 



Ca + 



9^ V3; 27 



Cf 



a + 0(a^ 



(B.3) 



The numerical evaluation of these as well as all the two loop corrections were given in |38] 
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C Expansions for various master integrals. 



In this section we present the expansions for several master integrals in powers of These 
were established by the methods of [31] which were discussed in section 3. First, we have for 
one mass scale 



i(i-/c2r, 0,0) 



9 3 



1 

7' 



IT I 

\- 1 In 

2 2 



3 1 , 

7 + 

4 4 



7^^ 



(C.l) 



For two non-zero masses we have 



i{WCAr,wcAr,o) 



vr + 2i — z In 



1 3, 
- + - In 


'Ca7^' 






2 2 


[ f^' \ 







(C.2) 



and 



/(!\/Cl7^-'\/Cl7',0) 



9 3^ V3 



+ 



2 



(C.3) 



which is real as expected. In the pure Gribov-Zwanziger case the main cases with three non-zero 
entries are 



i{iVCAr,wcAr,wcAr) 



and 



i{iVCAr,wcAr,-wcAY) 



47r 



9 3^ 



3 

Sir 3i , 

h 3i In 

2 2 



1 

1^ 



Cai' 



3 15, 

7 — rl'^ 

4 4 



Ca^ 



ISvri 



CaI^ 



6 



+ (C.4) 



9 3^ V3 



37r i 

\- 1 In 

2 2 



+ 



1 



In 



1 

Cat' 

13 57ri 
T ~ ~ 



+ 



(C.5) 



For either the Q or 7^ solutions we have 



1^ - f i 

9 3 V3 



+ 



A, 


In 




- 1 




In 




- 1 


+ 4 


In 


Mi 


- 1 
















[4 J 
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+ 4 



In 



lln 

2 



/4+ 



+ 



3 

+ 4 



In 



In 



+ 



4 
+ In 



I 2 2 



, 2 2 



In 


^1 


+ In 




- 1 


In 




+ ln 




- 1 



+ 



T. 



(C.6) 



for the case of three distinct non-zero masses. In each case the order symbols are intended to 
reflect the power of /i and the numerator factor therein is merely to have the correct overall 
dimensionful dependence. We have also given the expansion out to powers beyond that which 
we have indicated we are interested in for the overall vertex functions. This is because in the 
rearrangement of the numerator scalar products in the original integrals one can be left with 
terms such as p'^ and pq which are proportional to Hence, terms beyond the dimension 
four ones we are interested in for the overall vertex function need to be retained in the above 
examples and the expansion of the other basic master integrals. 
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